The Importance of Livestock Demography and Infrastructure in Driving Foot and Mouth Disease Dynamics.

Transboundary animal diseases, such as foot and mouth disease (FMD) pose a significant and ongoing threat to global food security. Such diseases can produce large, spatially complex outbreaks. Mathematical models are often used to understand the spatio-temporal dynamics and create response plans for possible disease introductions. Model assumptions regarding transmission behavior of premises and movement patterns of livestock directly impact our understanding of the ecological drivers of outbreaks and how to best control them. Here, we investigate the impact that these assumptions have on model predictions of FMD outbreaks in the U.S. using models of livestock shipment networks and disease spread. We explore the impact of changing assumptions about premises transmission behavior, both by including within-herd dynamics, and by accounting for premises type and increasing the accuracy of shipment predictions. We find that the impact these assumptions have on outbreak predictions is less than the impact of the underlying livestock demography, but that they are important for investigating some response objectives, such as the impact on trade. These results suggest that demography is a key ecological driver of outbreaks and is critical for making robust predictions but that understanding management objectives is also important when making choices about model assumptions.

Modeling nation-wide U.S. swine movement networks at the resolution of the individual premises.

The spread of infectious livestock diseases is a major cause for concern in modern agricultural systems. In the dynamics of the transmission of such diseases, movements of livestock between herds play an important role. When constructing mathematical models used for activities such as forecasting epidemic development, evaluating mitigation strategies, or determining important targets for disease surveillance, including between-premises shipments is often a necessity. In the United States (U.S.), livestock shipment data is not routinely collected, and when it is, it is not readily available and mostly concerned with between-state shipments. To bridge this gap in knowledge and provide insight into the complete livestock shipment network structure, we have developed the U.S. Animal Movement Model (USAMM). Previously, USAMM has only existed for cattle shipments, but here we present a version for domestic swine. This new version of USAMM consists of a Bayesian model fit to premises demography, county-level livestock industry variables, and two limited data sets of between-state swine movements. The model scales up the data to simulate nation-wide networks of both within- and between-state shipments at the level of individual premises. Here we describe this shipment model in detail and subsequently explore its usefulness with a rudimentary predictive model of the prevalence of porcine epidemic diarrhea virus (PEDv) across the U.S. Additionally, in order to promote further research on livestock disease and other topics involving the movements of swine in the U.S., we also make 250 synthetic premises-level swine shipment networks with complete coverage of the entire conterminous U.S. freely available to the research community as a useful surrogate for the absent shipment data.

Assessing intrastate shipments from interstate data and expert opinion.

Live animal shipments are a potential route for transmitting animal diseases between holdings and are crucial when modelling spread of infectious diseases. Yet, complete contact networks are not available in all countries, including the USA. Here, we considered a 10% sample of Interstate Certificate of Veterinary Inspections from 1 year (2009). We focused on distance dependence in contacts and investigated how different functional forms affect estimates of unobserved intrastate shipments. To further enhance our predictions, we included responses from an expert elicitation survey about the proportion of shipments moving intrastate. We used hierarchical Bayesian modelling to estimate parameters describing the kernel and effects of expert data. We considered three functional forms of spatial kernels and the inclusion or exclusion of expert data. The resulting six models were ranked by widely applicable information criterion (WAIC) and deviance information criterion (DIC) and evaluated through within- and out-of-sample validation. We showed that predictions of intrastate shipments were mildly influenced by the functional form of the spatial kernel but kernel shapes that permitted a fat tail at large distances while maintaining a plateau-shaped behaviour at short distances better were preferred. Furthermore, our study showed that expert data may not guarantee enhanced predictions when expert estimates are disparate.

Realistic assumptions about spatial locations and clustering of premises matter for models of foot-and-mouth disease spread in the United States.

Spatially explicit livestock disease models require demographic data for individual farms or premises. In the U.S., demographic data are only available aggregated at county or coarser scales, so disease models must rely on assumptions about how individual premises are distributed within counties. Here, we addressed the importance of realistic assumptions for this purpose. We compared modeling of foot and mouth disease (FMD) outbreaks using simple randomization of locations to premises configurations predicted by the Farm Location and Agricultural Production Simulator (FLAPS), which infers location based on features such as topography, land-cover, climate, and roads. We focused on three premises-level Susceptible-Exposed-Infectious-Removed models available from the literature, all using the same kernel approach but with different parameterizations and functional forms. By computing the basic reproductive number of the infection (R0) for both FLAPS and randomized configurations, we investigated how spatial locations and clustering of premises affects outbreak predictions. Further, we performed stochastic simulations to evaluate if identified differences were consistent for later stages of an outbreak. Using Ripley's K to quantify clustering, we found that FLAPS configurations were substantially more clustered at the scales relevant for the implemented models, leading to a higher frequency of nearby premises compared to randomized configurations. As a result, R0 was typically higher in FLAPS configurations, and the simulation study corroborated the pattern for later stages of outbreaks. Further, both R0 and simulations exhibited substantial spatial heterogeneity in terms of differences between configurations. Thus, using realistic assumptions when de-aggregating locations based on available data can have a pronounced effect on epidemiological predictions, affecting if, where, and to what extent FMD may invade the population. We conclude that methods such as FLAPS should be preferred over randomization approaches.

Effects of regional differences and demography in modelling foot-and-mouth disease in cattle at the national scale.

Foot-and-mouth disease (FMD) is a fast-spreading viral infection that can produce large and costly outbreaks in livestock populations. Transmission occurs at multiple spatial scales, as can the actions used to control outbreaks. The US cattle industry is spatially expansive, with heterogeneous distributions of animals and infrastructure. We have developed a model that incorporates the effects of scale for both disease transmission and control actions, applied here in simulating FMD outbreaks in US cattle. We simulated infection initiating in each of the 3049 counties in the contiguous US, 100 times per county. When initial infection was located in specific regions, large outbreaks were more likely to occur, driven by infrastructure and other demographic attributes such as premises clustering and number of cattle on premises. Sensitivity analyses suggest these attributes had more impact on outbreak metrics than the ranges of estimated disease parameter values. Additionally, although shipping accounted for a small percentage of overall transmission, areas receiving the most animal shipments tended to have other attributes that increase the probability of large outbreaks. The importance of including spatial and demographic heterogeneity in modelling outbreak trajectories and control actions is illustrated by specific regions consistently producing larger outbreaks than others.

Need for speed: An optimized gridding approach for spatially explicit disease simulations.

Numerical models for simulating outbreaks of infectious diseases are powerful tools for informing surveillance and control strategy decisions. However, large-scale spatially explicit models can be limited by the amount of computational resources they require, which poses a problem when multiple scenarios need to be explored to provide policy recommendations. We introduce an easily implemented method that can reduce computation time in a standard Susceptible-Exposed-Infectious-Removed (SEIR) model without introducing any further approximations or truncations. It is based on a hierarchical infection process that operates on entire groups of spatially related nodes (cells in a grid) in order to efficiently filter out large volumes of susceptible nodes that would otherwise have required expensive calculations. After the filtering of the cells, only a subset of the nodes that were originally at risk are then evaluated for actual infection. The increase in efficiency is sensitive to the exact configuration of the grid, and we describe a simple method to find an estimate of the optimal configuration of a given landscape as well as a method to partition the landscape into a grid configuration. To investigate its efficiency, we compare the introduced methods to other algorithms and evaluate computation time, focusing on simulated outbreaks of foot-and-mouth disease (FMD) on the farm population of the USA, the UK and Sweden, as well as on three randomly generated populations with varying degree of clustering. The introduced method provided up to 500 times faster calculations than pairwise computation, and consistently performed as well or better than other available methods. This enables large scale, spatially explicit simulations such as for the entire continental USA without sacrificing realism or predictive power.

Route optimization as an instrument to improve animal welfare and economics in pre-slaughter logistics.

Each year, more than three million animals are transported from farms to abattoirs in Sweden. Animal transport is related to economic and environmental costs and a negative impact on animal welfare. Time and the number of pick-up stops between farms and abattoirs are two key parameters for animal welfare. Both are highly dependent on efficient and qualitative transportation planning, which may be difficult if done manually. We have examined the benefits of using route optimization in cattle transportation planning. To simulate the effects of various planning time windows and transportation time regulations and number of pick-up stops along each route, we have used data that represent one year of cattle transport. Our optimization model is a development of a model used in forestry transport that solves a general pick-up and delivery vehicle routing problem. The objective is to minimize transportation costs. We have shown that the length of the planning time window has a significant impact on the animal transport time, the total driving time and the total distance driven; these parameters that will not only affect animal welfare but also affect the economy and environment in the pre-slaughter logistic chain. In addition, we have shown that changes in animal transportation regulations, such as minimizing the number of allowed pick-up stops on each route or minimizing animal transportation time, will have positive effects on animal welfare measured in transportation hours and number of pick-up stops. However, this leads to an increase in working time and driven distances, leading to higher transportation costs for the transport and negative environmental impact.

High frequency of functional extinctions in ecological networks.

Intensified exploitation of natural populations and habitats has led to increased mortality rates and decreased abundances of many species. There is a growing concern that this might cause critical abundance thresholds of species to be crossed, with extinction cascades and state shifts in ecosystems as a consequence. When increased mortality rate and decreased abundance of a given species lead to extinction of other species, this species can be characterized as functionally extinct even though it still exists. Although such functional extinctions have been observed in some ecosystems, their frequency is largely unknown. Here we use a new modelling approach to explore the frequency and pattern of functional extinctions in ecological networks. Specifically, we analytically derive critical abundance thresholds of species by increasing their mortality rates until an extinction occurs in the network. Applying this approach on natural and theoretical food webs, we show that the species most likely to go extinct first is not the one whose mortality rate is increased but instead another species. Indeed, up to 80% of all first extinctions are of another species, suggesting that a species' ecological functionality is often lost before its own existence is threatened. Furthermore, we find that large-bodied species at the top of the food chains can only be exposed to small increases in mortality rate and small decreases in abundance before going functionally extinct compared to small-bodied species lower in the food chains. These results illustrate the potential importance of functional extinctions in ecological networks and lend strong support to arguments advocating a more community-oriented approach in conservation biology, with target levels for populations based on ecological functionality rather than on mere persistence.